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ABSTRACT 

The faint stellar halos of galaxies contain key information about the oldest stars and 
the process of galaxy formation. A previous study of stacked SDSS images of disk 
galaxies has revealed a halo with an abnormally red r — i colour, seemingly inconsistent 
with our current understanding of the stellar populations inhabiting stellar halos. 
Measurements of this type are however plagued by large uncertainties which calls for 
follow-up studies. Here, we investigate the statistical properties of the faint envelopes 
of low surface brightness disk galaxies to look for further support for a red excess. A 
total of 1510 nearly edge-on, bulgeless low surface brightness galaxies were selected 
from the SDSS Data Release 5, rescaled to the same apparent size, aligned and stacked. 
This procedure allows us to reach a surface brightness of \i r ~ 31 mag arcsec -2 . After 
a careful assessment of instrumental light scattering effects in the stacked images, we 
derive median and average radial surface brightness and colour profiles in g,r and i. 
The sample is then divided into 3 subsamples according to g — r colour. All three 
samples exhibit a red colour excess in r — i in the thick disk/halo region. The halo 
colours of the full sample, g — r = 0.60 ± 0.15 and r — i — 0.80 ± 0.15, are found to 
be incompatible with the colours of any normal type of stellar population. The fact 
that no similar colour anomaly is seen at comparable surface brightness levels along 
the disk rules out a sky subtraction residual as the source of the extreme colours. A 
number of possible explanations for these abnormally red halos are discussed. We find 
that two different scenarios - dust extinction of extragalactic background light and 
a stellar population with a very bottom-heavy initial mass function - appear to be 
broadly consistent with our observations and with similar red excesses reported in the 
halos of other types of galaxies. 

Key words: galaxies: halos - galaxies: stellar content - galaxies: structure - galax- 
ies:photometry - stars: low-mass - diffuse radiation 



1 INTRODUCTION 

The stellar halos of galaxies are nearly spherical, old and 
metal-poor structures with little net rotation. Due to the 
long orbital time scales of halo stars, these halos contain the 
fossil record of galaxy assembly, and observations of resolved 
halo stars in the Milky Way, Andromeda and other nearby 
galaxies (e.g. Mouhcine et al.|2005 Ibata et al.|2007 Helmi 



2008} |Rejkuba et al.|2009||de Jong et al.|2009| > have provided 
a wealth of information on such structures. Current simula- 
tions of st ellar halos (e.g. [Bullock fc Johnston|[2005| |Helmi| 
et al.|2006[|Gauthier et al.|2006||Purcell et al.|2007||De Lucia| 



fc Helmi|2008| |Font et al.|2008|jZolotov et al.|2009| |Stmson] 



et al. 2009) indicate that they form through a combination 
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of disruption of accreted dwarf galaxies, and by migration of 
old stars outwards from the inner regions of galaxies. Mod- 
els of this type are able to explain many of the features 
observed - including the numerous streams from disrupted 
satellite galaxies, the inner-outer halo dichotomy, and the 
existence of stellar halos around dwarf galaxies. However, 
there is a lingering uncertainty concerning the integrated 
colours of halos. 

Observationally, there are two complementary tech- 
niques for studying the halos of galaxies: direct star counts 
and surface photometry. Star counts can trace the structure 
of stellar halos out to much larger distances (or, equiva- 
lently, to fainter surface brightness limits) than surface pho- 
tometry can, but is currently only applicable for galaxies 
at distances out to ~ 10 Mpc. Surface photometry (i.e. the 
study of the integrated light from large numbers of unre- 
solved sources within the each system) can be used for far 
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Figure 1. Characteristic spectra of the galaxies in the three subsamples A,B and C obtained from the SDSS. 



more distant galaxies and in principle gives a more direct 
measurement of surface brightness, since it is sensitive both 
to stars at masses below the detection threshold of current 
star counts and to diffuse light from the interstellar medium. 
This is, on the other hand, a far more challenging technique, 
which requires that the sky flux can be subtracted with ex- 
cruciating accuracy and that various instrumental scatter- 
ing effects can be kept under control. Curiously, a number 
of attempts to study the halos of galaxies through surface 
photometry have revealed halo colours that are too red to 
be straightforwardly reconciled with the halo populations 



captured through the study of resolved stars (see Zackrisson 



|et al.|2006| and|!Zackrissoii ct al. 2007]for a general discussion 
of this phenomenon). 

The first indication of a red excess in the halos of disk 
galaxies was found through BVR1JK photometry of the 
edge-on system NGC 5907. The interest was initiated by 
the discovery of a faint extended i?-band halo around this 
galaxy ( Sackett et al.|1994 1 , which triggered follow-up stud- 



ies in the optical and near-IR (Lequeux et al. 1996 Rudy 



et al. 1997 James & Casali 19981. A common interpreta- 



tion of the red excess was that the galaxy halo hosted an 
unusually large population of low-mass stars. Later, deep 



images in the optical (Shang et al. 1998 Zheng et al. 1999 



|Martmez-Delgado et al.||2008[ ) revealed complex arc struc- 
tures in the halo of the galaxy that, when mixed with scat- 
tered light from stars in the field, may have caused confusion 
in the analysis of previous data of lower S/N. Therefore the 
reality of a true, diffuse halo in NGC 5907 at the surface 
brightness levels where the red excess had been measured 
was questioned by |Zheng et al.| (|1999||. St ill, an HST study 
of resolved stars in the halo (Zcpf et al. 20001 indicated a 



relative paucity of massive stars which could be interpreted 
either as a very low metallicity of the stellar population or 
as a bottom heavy initial mass function (IMF), inconsistent 
with the halo colours as observed from the ground. 

Surface photometry at wavelengths outside the optical 
at low brightness levels are valuable as a complement to the 
optical. If the red excess is of thermal origin, it is likely to 
show up also when compari ng the optical to the nea r-IR 
fluxes. A few years ago (see Bergvall & Ostlin 20021, we 



obtained deep optical/near-IR images of a handful of Blue 
Compact Galaxies (BCGs). When inspecting the faint halo 
regions we found a strong red excess in the near-IR, again 
inconsistent with normal stellar populations. Dust emission 
or extinction could be ruled out as a major cause of the ex- 
cess in our sample. The infrared excess was later confirmed 



for a larger sample of BCGs ( Bergvall et al.|[2005 1 . The re- 
sults from these studies relied on contemporary models of 
galaxy spectral evolution. An uncertain ingredient in these 
models is the late stages in the evolution of late type stars. 
Recent models predict redder colours of stars of intermedi- 
ate masses which tends to reduce the discrepancy between 
observations and models (Zackrisson et al. 20091. Still, we 



are left with an uncomfortable fine tuning of age of the halo 
stars and either a metallicity or dust reddening that must 
be higher than in the starburst region. 

One problem with all these investigations is the diffi- 
culties of dealing with data at very low surface brightness 
levels, much below 1% of the background sky. Under such 
conditions it becomes critical how the subtraction of the 
emission from the sky and scattered light is performed. This 
is a problem that we can never overcome as long as we use 
data from groundbased instruments. A second problem is 
due to the low photon flux rate in the studied regions, caus- 
ing an uncomfortably large photometric inaccuracy. If we 
want to study the halos of individual galaxies, the only way 
we can avoid this problem is to use very large telescopes or 
extend the integration times to the extreme. 

An alternative approach was tested by |Zibetti et al.| 
(20041, hereafter Z04. It is based on the statistical analy- 



sis of a large sample of galaxies with common fundamental 
properties (essentially morphologies, luminosities and stel- 
lar content). In such a fairly homogeneous sample one may 
assume (since one cannot measure each individual galaxy 
with sufficient S/N) that the physical properties of the faint 
halos of the galaxies are also similar. An increase in S/N can 
therefore be obtained by stacking all images after rescaling 
and aligning the individual images of the galaxies. Based 
on data from the Sloan Digital Sky Survey (SDSS), Z04 
tested this approach and digitally stacked images of 1047 
nearly edge-on disk galaxies. As a result they managed to 
reach extremely low surface brightness levels, fi r ~ 31 mag 
arcsec -2 . When measuring the colours in the polar direc- 
tion, they found that the colours (g — r= 0.65±0.1, r — i— 
0.60±0.1) of the faintest regions where reliable data could be 
obtained, were deviating in r — i with 2a or about 0.2 m from 
the reddest population observed in areas of higher surface 
brightness. Although acknowledging that a small amount of 
the red excess could be credited to light scattering effects, 
they still concluded that "r — i [...] is difficult to reconcile 
with any theoretical models, even allowing for ad hoc mod- 
ified IMFs dominated by low-mass stars and high metallic- 
ity" . One might suspect that the excess they found in the 
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optical (i band) was a confirmation of the red excess we 
found in the near-IR. This remains to be confirmed, how- 
ever (see discussion in e.g. |Zackrisson et aL 20091. While 



Table 1. Properties of subsamples A, B, C and the total sample, 



being confident with the existence of an abnormal red ex- 
cess, Z04 point out that their "colour measurements provide 
inconclusive and troublesome results". Recently, there was 



a claim by de Jong ( 2008 1 that the abnormal colours could 
have been caused by an improper treatment of the point 
spread function (PSF) when analysing the effects of light 
scattering. An independent confirmation or falsification of 
the results by Z04 is therefore obviously of highest concern. 

In the present investigation we carry out a follow-up of 
the work of Z04. Preliminary results were presented by |Cald-] 
well & Bergvall ( 2007 1. The sample we select for our study is 



different from that of Z04 and includes more galaxies, thus 
increasing the S/N after stacking. It focuses on another type 
of galaxy, Low Surface Brightness Galaxies (LSBGs), and 
thus allows us to make more general conclusions about the 
phenomenon, if a red halo is found. We chose to work with 
LSBGs mainly for two reasons. Firstly, they have a low dust 
content (e.g Bergva fl et al.|1999 Holwerda et al.|2005 Rah 



|man et al.|2007[ ), thereby essentially removing the potential 
problem of dust reddening. Secondly, the influence of the 
light scatter from the bulge is less severe than in normal 
disks. This second argument is weaker however, since we do 
not know a priori if the red halo perhaps is physically cor- 
related with the bulge component. An additional advantage 
is that we could learn about the properties of the spheroidal 
component of LSBGs, hitherto almost unknown. 

The method we will use is similar to that used by Z04 
with a few modifications as described below. We select im- 
ages of LSBGs from SDSS in the g, r and i bands, scale and 
align and finally stack the images. In the following we will 
describe the selection criteria and reduction procedures in 
more detail. 



2 DATA ANALYSIS 
2.1 Selection criteria 

We selected a sample of nearly edge-on LSBGs from the 



SDSS Data Release 5 (DR5) database ( Adelman-McCarthy 
12007) . 

In the selection of targets we used the following criteria: 

(i) 23< jU r ,50 <25 mag arcsec -2 , where /v,so is the mean 
surface brightness inside the radius encompassing 50% of the 
light in the r-band 

(ii) i-band isophotal diameter of a ^ 10 arcsec 

(iii) i-band axial ratio of b/a < 0.25 

(iv) redshift z < 0.2 

The red i band was chosen in order to minimize the 
effects on the size parameters due to patchiness caused by 
young blue star forming regions and dust extinction. The 
axial ratio was chosen in order to avoid too much contribu- 
tion from the inclined disk to the halo light. We made the 
choice with the aim to optimize the S/N of the halo light. 
The halo region of interest that we will discuss below starts 
at about 60 pixels in the vertical direction. At that point 
the inclined disc contributes about 20-30% to the light in 
the i band. At larger distances the relative contribution is 



Name 


No. of galaxies 


Colour constraint 


A 


447 


g — r < 0.55 


B 


655 


0.55 < g-r < 0.75 


C 


408 


0.75 < g-r 


T 


1510 





considerably lower. The useful region for determination of 
the halo colour with the present axial ratio limit is ~ 80-110 
pixels. Should we increase the limit, we would obtain more 
galaxies but the usefiul region would shrink and the reverse 
would be the case if we should decrease the limit. In the end 
we are limited by the value of the intrinsic flattening. 

The final sample contained 1510 galaxies. We used Pet- 
rosian magnitudes and divided the sample by g — r colour 
into three subsamples, a "blue" sample (A), a "green" sam- 
ple (B) and a "red" sample (C). Set Table 1 for details. 
The goal of this subdivision was to search for possible cor- 
relations between the strength of the red excess and the 
total colour /luminosity. The maximum redshift was chosen 
to avoid uncertain fc-corrections, uncomfortably noisy data 
and to permit a direct comparison with the study by Z04. 
In fact very few galaxies fulfilling the remaining criteria are 
found at higher redshifts. 

All images were corrected for atmospheric extinction 
based on primary standards observed regularily at Apache 
Point (where the SDSS observations are carried out). Galac- 
tic extinction was corrected for according to Schlcgcl et al. 
(1998). No fc-corrections were applied. This will affect the 
red sample the most since those galaxies have the largest 
median redshift. The correction is however relatively small 
for a normal galaxy of intermediate/late type. At the me- 
dian redshift (< z >~ 0.1) of the reddest sample, the k- 
correction is <0.2 m in both colours (see e.g. Fukugita et al. 



19951. A fc-correction of the halo colours is anyhow useless 
because, when it comes to regions showing a red excess, we 
do not know a prion the origin of the radiation. In section 
4 we will discuss how the colours are affected at different 
redshifts assuming various types of radiation sources. 

2.2 Sample properties 

The following figures show comparisons between the proper- 
ties of each subsample and the full sample. Typical examples 
of spectra of galaxies from the three subsamples are shown 
in Fig. [I] They show that there is a large difference in age of 
the stellar population dominating the light within the SDSS 
aperture (0=3") from a few 10 7 years to several Gyr. 

The distribution of redshifts is shown in Fig. [2] There 
is a clear trend towards higher redshifts as we go from the 
blue to the red sample, with a median redshift of z~0.06. 
The correlation between redshift and colour reflects the fact 
that the galaxies in the blue sample are smaller and there- 
fore reach the limiting apparent diameter at smaller red- 
shifts than the redder ones. There is also a related correla- 
tion between colour and absolute luminosity shown in Fig. 
[3] displaying the distribution of absolute magnitudes. The 
bluer Sample A falls within the classical range of typical LSB 
galaxy luminosities of -13 ^ Mb ^-19 ( Impey et al. 1(1996 1. 
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Figure 2. The redshift distribution of the subsamples and the 
full sample. 




Figure 3. The distribution of absolute magnitudes in the i band 
of the galaxies in the subsamples and the full sample. 
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Figure 4. The distribution of exponential scalelengths of the 
galaxies in the subsamples and the full sample. 
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Figure 5. The two-colour distribution of the galaxies in the three 
subsamples, separated at colours g — r=0.55 and 0.75. 



These galaxies are normally gas rich and globally star form- 
ing. The peaks of the green and red samples are moving more 
toward normal spiral galaxy type absolute magnitudes. 

Fig. [4] presents the exponential scalelength distribu- 
tions. The range of scalelengths is from 1 to 16 kpc h~ 
with a median value of 5.5 kpc h . This is somewhat larger 
than that of high surface brightness disk galaxies studied by 
|Dalcanton et al.| |T997| ; |MacArthur et al.] ( |2003[ > and twice 
the size of the Z04 sample. It confirms and strengthens the 
results from a previous study of LSB galaxies by |Zwaan et al.| 
( |1995[ ). In this study, the difference in scale length was in- 
terpreted as consequence of a larger spin for a given mass in 
the LSB sample. 

The colour-colour distribution is shown in Fig. [5] The 
scatter of the whole sample is most apparent in the blue 
subsample and most constrained in the redder LSB galaxies. 
This may be explained as an effect of the patchiness of the 
blue LSB galaxies due to a larger fraction of star forming 
regions. 



2.3 Stacking preparations 

Since our main result, the red excess, could be regarded as 
controversial, we want to make sure that our data retrieval 
and processing can be reproduced by anyone who reads this 
article. In the following we will therefore describe the pro- 
cedure in detail. 

The galaxy target field images were obtained from the 
SDSS archive. The resulting raw SDSS images are of size 
2048x1489 pixel 2 or 13.5x9.8 arcmin 2 . Each of these frames 
were bias subtracted and flat field corrected. For each galaxy 
we obtained the galaxy centres, position angles, scalelengths 
and atmospheric extinction corrected magnitudes from the 
SDSS database. The so prepared images were then passed 
through a number of processes one by one where they were 
masked, shifted and rescaled before finally being stacked, 
averaged and median filtered, all with the combined use of 
the ESO MIDAS command language, version 08FEB pll.l on 
MAC/OSX and S Extractor v2.4.4 (Source Extractor, [Bertin| 
fc Arnouts|1996 1 packages in the following way. 
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First SExtractor was used to automatically remove all 
objects in the frame except for the target galaxy itself. A 
first fit to the sky background was initially made based on 
a course grid of resolution 512x512 pixels. Since the typ- 
ical diameter of a target galaxy is about 10% of the grid 
size we did not expect serious problems with fluctuations 
in the fitted sky on the size of the halo we were exploring. 
After the fit, the background was subtracted and all ob- 
jects, defined as at least 5 connected pixies with fluxes ^ 
1.5 times the local pixel mean error in the sky level, were 
automatically analysed, identified as stars or galaxies and all 
objects except for the target galaxy were then flagged into 
a separate file. The galaxy was now the only object (de- 
fined as above) present in the frame. The resulting dummy 
file had the galaxy moved to its center and rotated so that 
the major axis of the galaxy became horizontal. The image 
was then scaled to a standard exponential scalelength of 20 
arcseconds while conserving surface brightness. The central 
1024x1024 pixels were extracted into another file: the im- 
age subfile. Because the first sky fit was made with a coarse 
grid over a region many times larger than the galaxy, weak 
local gradients at scales below the size of the mesh could 
still be present. Therefore a second sky subtraction was per- 
formed in the smaller field, based on a grid with a mesh size 
of 256x256 pixels, excluding the flagged regions and the tar- 
get. The mesh size was thus half the size used for the first fit 
to the sky background. New objects were now identified and 
flagged. To make sure that the wings of the smaller stars and 
galaxies were completely removed we flagged 5 extra pixies 
around the flagged regions around all detected objects. 

All further pruning was done by eye to identify objects 
not suitable for stacking. Some false detections were obvi- 
ous as problems with the SDSS PHOTO reduction pipeline. 
Stars superposed on the stacked galaxy that were not iden- 
tified by SExtractor were flagged. This process is somewhat 
subjective since at low relative brightnesses it is difficult to 
distinguish between a superposed star in the Milky Way and 
a luminous star cluster in the galaxy. Potential target galax- 
ies were rejected if they had obvious warps, excessive dust 
lanes, bright bulges, were in interaction with nearby com- 
panions or showed other irregular morphologies. A further 
set was discarded due to the closeness of nearby bright stars. 
Fainter structures were identified after filtering the images 
with a gaussian filter of a — 5 pixels. The filtering makes 
it easier to detect wings from bright stars just outside the 
field borders and other low frequency variations. Unwanted 
regions were thus flagged or the whole image was deleted. 
The returned rough sample was culled to a total of 1510 
galaxies: 447 in sample A, 655 in sample B and 408 in sam- 
ple C. The relative fraction of galaxies fulfilling the selection 
criteria that were finally included in the stacking were 62%, 
71% and 69% for the A,B and C samples respectively. 

2.4 The image stacking 

The technique we use for image stacking largely follows that 
used by Z04. The stacking process aims to increase the 
signal/noise (S/N) by median filtering a large number of 
rescaled and aligned images of assumingly conforming galax- 
ies. The difference between our procedure and Z04, is that 
we use both averaging and median filtering while they use 
the mode, based on an empirical relation. If there is a wide 
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Figure 6. The PSF profile in g (dashed), r (dotted) and i (solid) 

diversity in morphologies and luminosity gradients one may 
worry about what the result will tell us, at least if we use me- 
dian filtering. Therefore it is important to make comparisons 
between the results from median filtering and averaging. We 
will show below that from this aspect, the result we obtain 
is quite robust. 

The image stacking technique consists of three basic 
steps: (1) masking out all unwanted surrounding objects (2) 
aligning, rotating and rescaling the LSBG images to make 
them superposable (3) combining the images to produce an 
average and median. All the above steps were combined 
in an in-house developed programme written in MIDAS. 
This programme processes the SDSS images through the use 
of the MIDAS and SExtractor software packages and local 
FORTRAN routines. The median and average image of the 
resulting individual masked, centered and rescaled images 
were then computed using MIDAS routines. The S/N within 
the areas we study of each individual galaxy is proportional 
to the square root of the number of photons received, i.e. 
the surface area. In the stacking procedure, the images were 
therefore given weights proportional to the original scale- 
size squared. This has the consequence that the integrated 
weight given galaxies at a certain distance is independent 
of distance. Thus, nearby galaxies are strongly overrepre- 
sented in the stacked image. The ratio between the largest 
and smallest galaxies in the subsamples is typically ~ 10. 
However, no single large galaxy dominates any of the three 
subsamples. After the stacking, a flat surface fit to the back- 
ground sky in 15 boxes around the galaxy was made and 
subtracted. A final background correction was carried out 
in connection to the integration of the luminosity profiles 
(see below) . 

2.5 The error in the chosen sky level 

As a check to find out whether the sky noise is purely Pois- 
sonian or structured on some specific scale we carried out a 
simple test. We compared the pixel-to-pixel standard devia- 
tion Usky calculated from data obtained in boxes of various 
sizes, up to approximately the scale of the image, in our case 
avoiding the corners that contain more stellar residuals than 
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Figure 7. Isophotal contour plots of the stacked images of the subsamples (A, B and C) and the total sample (T). The images have 
been smoothed with a gaussian filter with <r=7 pixels. The contour levels arc plotted in steps of 1 magnitude. The faintest level in the 
g, r and i images of each sample is A) 30.1, 29.8, 30.1 B) 30.0, 29.7, 29.2 C) 29.5, 29.2, 28.7 T) 30.6, 30.1, 29.7 mag arcscc~ 2 respectively. 
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the rest of the image. In a completely structureless (white 
noise) sky, a s ky should be independent of the size of the box. 
We found this to be close to fulfilled so apparently the me- 
dian filtering and the simple sky subtraction was sufficient 
to remove measurable correlated structures outside the im- 
age of the galaxies. The error in the background zeropoint 
level (after subtracting the sky values from the data) was 
therefore set equal to the mean error of the mean of the 
median values of the pixels in the 15 boxes used for the sky 
subtraction. 



2.6 The Point Spread Function 

A serious problem when dealing with weak signals in images 
is the influence of scattered light from the target and objects 
in the environment. The importance of this potential prob- 
lem can be estimated if the shape of the PSF is known. This 
problem was discussed at length in the investigation by Z04. 
To tackle the problem, the authors modelled the PSF in the 
different filters by a combination of a gaussian core plus an 
exponential wing. They concluded that the influence of the 
PSF on the data in the red halo region was insignificant. 

Here we have chosen a more direct approach. To ob- 
tain the PSF, we stacked bright stars, one from each galaxy 
field after cleaning and rejecting bad images in the same 
manner as before. The result in shown in Fig. [6] Obviously 
there is a clear difference in the PSF between the i window 
and the g and r windows. The PSFs of SDSS images have 
independently been measured by |Wu et ak] ( |2005[ ) and [de] 
Jong (20081. Our results agree very well with theirs down 
to the levels that might be relevant in this investigation, i.e. 



about 15 magnitudes below peak surface brightness. The 
PSF profiles have been derived with the same procedure as 
the galaxy profiles, and the tight agreement between our 
results and those of Wu et al. and de Jong supports the ro- 
bustness of our method. The shape of the PSF is markedly 
different from that derived by Z04. The one we derive is sig- 
nificantly more extended and clearly differs between i and 
the other filters at distance from the centre that in the mean 
is close to the distance at which Z04 find the red colour ex- 
cess. This potential problem is discussed in detail in the 
paper by de Jong. 



2.7 Correcting for the PSF 

2.7.1 Deconvolution 

We used two independent methods to correct for the de- 
grading of the images due to light scattering, as quantified 
by the PSF. The first one was simply to deconvolve each 
individual image before stacking. We carried out a decon- 
volution of the full sample using the MIDAS software. The 
deconvolution was done before the rescaling and stacking of 
each individual image using the task DECONVOLVE in MI- 
DAS, based on the iterative algorithm by L.B. |Lucy| ( |1974[ ). 
Before the deconvolution, all images were rebinned to twice 
the original stepsize in order to avoid long processing times. 
The size of the PSF image after rebinning was 95x95 pix- 
els. The deconvolution increases the high frequency noise 
considerably due to the influence of the background objects 
and may cause unwanted ringing effects in the data. The 
problems are most severe in the deconvolution of galaxies 



of small apparent sizes. When we inspected the results from 
the deconvolution we found structures at low fluxes that ap- 
peared artificial. We will therefore not discuss these results 
other than when it concerns data with weak radial gradients 
at flux levels safely above the background noise. 



2.7.2 Using a model of the luminosity distribution as a 
template 

The other method we used to correct for the PSF is based 
on a model of a typical sample galaxy. A similar procedure 
was recommended (although not practised) by de Jong as a 
proper way to account for the PSF problems in the Z04 sam- 
ple. Our model was constructed in the following way. The 
35 apparently largest galaxies were selected, deconvolved, 
scaled to the same apparent scale sizes and stacked. The de- 
convolution seemed to work well for these galaxies, having 
large apparent sizes. The resulting stacked image was used 
as a template to a model that was derived iteratively from a 
mixture of disks and flattened spheroidals with different Bat- 
tenings, radial surface brightness dependencies and relative 
peak fluxes. A good fit to the template galaxy was obtained 
from a simple mixture of a spheroidal with a flattening of 
b/a=l/3-l/2 and an r 1 / 4 surface brightness dependence su- 
perposed on an exponential disk. The exponential disk had 
two sections with different slopes, i.e. a "truncated disk"-like 
morphology. No bulge component was necessary to obtain a 
satisfactory agreement between the model and the template. 
The structural parameters of the model were adjusted to op- 
timize the fit to the template. In the end the deviations in 
surface flux in regions of the size of the seeing disk were 
less than about 15%. The parameters of the final model are 
shown in Table 2. For this subsample one pixel corresponds 
to ~ 0.32 arcsecs. This model was scaled to the same size 
as each galaxy in the sample and then treated equally in 
the stacking procedure. In a parallel sequence we convolved 
the model with the PSF after scaling and then stacked the 
images. In the end we thus obtained one image produced by 
the pure template images and the other by the PSF down- 
graded images. The ratio between these images could then 
directly be used to correct the images of the sample galaxies. 
The procedure seems to be quite independent of the model, 
except for the central region where the luminosity gradients 
are large. As an additional test we simply used the stacked 
image as a model, i.e. based on already degraded images. 
The result of the halo colours was very similar to the more 
appropriate method described above. The PSF corrections 
mainly affect the fainter parts of the r — i colour profile, 
making it bluer after the correction (see also Sect. 3.3). In 
the discussion below we will mainly refer to data corrected 
for the PSF using this method. 



3 RESULTS 

3.1 Luminosity distribution 

The stacked images are shown as contour plots in Fig. [7| A 
slight increase in the noise and the mean level of the sky is 
seen towards the edges. This is due to influence of objects 
outside the frame limits where the flagging procedure was 
not completely successful. The typical enhancement of the 
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Table 2. Model galaxy photometric parameters. 



Filter Inner disk 




Outer disk 




Truncation radius 


Flattened spheroidal 


Mo h 


b/a 


MO h 


b/a 


rtrunc 


MO r e c/ a 


mag arcsec -2 pixels 




mag arcsec -2 pixels 




pixels 


mag arcsec -2 pixels 


g 23.81 41 


0.12 


20.55 18 


0.12 


95 


21.92 132 0.34 


r 23.42 48 


0.16 


19.80 17 


0.16 


86 


21.13 146 0.34 


i 23.37 51 


0.16 


22.06 17 


0.16 


86 


20.39 100 0.50 



Mo ; Central surface brightness 
h: Exponential scalelength 

b/a: Observed (projected) minor/major isophotal axes ratio 
r e : Effective radius 



sky region due to this is ~10% of the signal at the faintest 
reliable part of the luminosity profile. As is seen in the plots, 
the background surface in the region close to the galaxy is 
essentially flat. The zero-point error of the sky background 
estimated from the statistics of the boxes used for the sky 
subtraction is typically about 1% of the surface brightness 
in the faintest regions where we consider the photometry to 
be reliable. 

We find no significant difference in the morphologies of 
the three subsamples from the contour maps. Small differ- 
ences between the samples can be found if we inspect the 
result we obtain by dividing the images in one sample with 
the corresponding images in another sample. We then find 
that the C sample has a more prominent contribution from 
the central region compared to samples A and B. 



3.2 Luminosity profiles 

We will now discuss the luminosity and colour profiles of the 
stacked images as derived from three regions of the images. 
We will call them halo, thick disk and thin disk. The halo is 
the faintest and most extended region of the stacked images. 
It is so faint that its structure is difficult to determine. Its 
presence is mainly revealed in the 1-D luminosity profiles. 
We will assume that it is spherical. The thick disk is defined 
by the faintest regular isophotes of the stacked images. It is 
close to ellipsoidal in shape but must be significantly flatter 
than the halo. There may be a smooth transition between 
the thick disk and the halo. The thin disk is the dominating 
morphological structure, defined by the faint isophotes of 
the individual images. It is clearly separated from the thick 
disk/halo in the r — i colour profiles. More details are given 
below. 

The halo plots were obtained from data in two opposite 
wedges along the polar axis of opening angle of 60°. The 
halo profile is measured in circular arc strips. These mea- 
surements correspond to those of Z04. The equatorial data 
along the thin disk were measured in elliptical strips in ho- 
risontally oriented wedges with an opening angle of 30° . The 
median inclination of a sample of infinitely thin disks with 
random orientations restricted to an observed minor/major 
axis ratio of b/a <0.25 (one of the selection criteria) is ~ 
81° or fe/a=0.16. The measured axis ratio of the i image 
of the full sample, based on the isophote encompassing the 
faintest regions of the thin disk, is b/a = 0.24, indicating a 
finite thickness of the disk. These data can be used to cal- 
culate the intrinsic minor/major axes ratio of the disk, p, 



from 



(b/a) 2 



where i is the inclination (Hubble 1926 Holmbcrg 1950J). We 
then obtain p ~ 0.18 . This is slightly flatter than what was 



obtained by Yoachim & Dalcanton ( 2006 1 for high surface 



brightness disk galaxies. We have used the observed axis 
ratio of 0.24 when we integrate the luminosity profile of the 
disk. 

As we discussed in the previous section, we obtained a 
satisfactory fit to the template galaxy using a model galaxy 
based on two components only - a truncated thin disk and 
a flattened spheroidal. At least at these isophotal levels 
(only 35 galaxies were used for the template), no bulge or 
halo component was needed. As we will show below, when 
we analyse the total sample (1510 galaxies) this conclusion 
largely persists. One must keep in mind however that our 
model of the structural properties is quite simple. In the 
faint region the stacked images seem to fit well to a flat- 
tened spheroidal but we cannot exclude that an equally good 
fit would be obtained with a mixture of a thick exponen- 
tial disk and a more extended spheroidal halo component. 
For the sake of simplicity, we call this dominating flattened 
structure the thick disk. It has a flattening in the PSF cor- 
rected image in the i band of b/a = 0.39±0.01, based on the 
three subsamples. The corresponding intrinsic flattening is 
0.36. Thus the ratio of the flattenings between the thick and 
thin disk is 2.05, slightly lower than the corresponding value 
for normal disks of 2.35 ( Yoachim fc Dalcanton|[2006[ ). LSB 
galaxies normally have flatter disks than normal galaxies 
which generally is interpreted as an evolutionary effect dur- 



ing the slow collapse of the disk (Kregel et al. 20051. The 



observed flattening of the thick disk is significantly smaller 
in the g and r bands, 0.26 and 0.29 respectively. This condi- 
tion makes the i band flux become more prominent relative 
to the other bands along the polar direction. We call this 
phenomenon the "red excess" in the discussion below. The 
luminosity/colour profiles of the thick disk (see below) are 
based on data collected in vertically oriented wedges with 
an opening angle of 120°. We integrated in elliptical strips 
corresponding to a flattening of 0.39. 

The luminosity profile of the total sample along the po- 
lar direction is shown in log-log and log-lin form in Fig. 
[8] The corresponding equatorial profiles along the disk are 
shown in log-lin form in Fig. [9] Errors are based on statistics 
of the sky counts in the strip used to derive the luminosity 
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Figure 8. Luminosity profiles along wedges of 60° opening angles in the polar direction of the stacked images of the total sample (T). 
The bars indicate mean errors as described in section 3.2. Profiles are presented in logarithmic and linear radial scale. 




Radius (pixels) 

Figure 9. Luminosity profiles along the major axis of the thin 
disk in sample T. The opening angle is 30°. A disk flattening of 
b/a=0.24 is assumed. The bars indicate mean errors. The physical 
scale in this and the following diagrams based on the full sample 
is approximately 500 pc pixel . 



profile at that specific radius. Thus the mean error in each 
bin was calculated from the mean error in the mean of the 
pixel values to which was added the estimated mean error 
in the sky level in quadrature. All data in the diagrams have 
been corrected for the effects of the PSF. 

The disk profile in i follows an exponential distribution 
in surface brightness. Subtracting this component would 
leave a structure resembling an extended bulge component. 
However, as we show in our model, it can be explained by 
the dominance of the thick disk in the central region. In the 
polar direction the curves display a behavior typical of a 
Sersic profile with n>l except in the outermost parts where 
they tend to be more exponential-like. 



3.3 Colour profiles 

We will now discuss to what extent we can find any support 
for a red colour excess at the so faint surface brightness 



levels. We will start by discussing the thin disk colours and 
then go to the thick disk and halo data. 



3.3.1 Thin disk 

Fig. [10] shows the equatorial colour profiles of the thin disk 
of the full sample. After correction for a small amount of 
dust reddening, the central colours are in agreement with a 
middle-aged/old stellar population. The corresponding age 
and metallicity depend on the amount of reddening assumed. 
The colours then drop smoothly to bluer colours and end up 
at g — r ~ 0.5, r — i ~ 0.2, corresponding to a normal stel- 



lar population with a standard IMF (cf. Fig. 151. An age- 
metallicity degeneracy prohibits us from being more specific 
about the age, but most likely the metallicity is low in the 
outer disk which would mean a high mean age of the popu- 
lation. We see no significant trends in the colours becoming 
redder at low surface brightness levels. 



3.3.2 Thick disk and halo 

We begin our discussions with the subsamples A, B and C. 
Fig [Tl] shows the colour profiles of the thick disk. Here the 
colours become somewhat bluer as we go from the centre 
outwards. At about 100 pixels they suddenly rise to a very 
red colour in r — i while they do not change significantly 
in g — r. The change occurs at a distance corresponding to 
the edge of the thin disk and is strongest in the B sample. 
There the colours become extreme, g — r — 0.4T0.10, r — i — 
1.2T0.15. It is exceedingly difficult to reconcile these colours 
with a normal stellar population, even if we allow the slope of 
the IMF to take extreme values. Compared to normal galax- 
ies in the SDSS, the disagreement is extreme. For instance, 
galaxies of the reddest morphological type in the SDSS (E 
galaxies) have g 



=0.83±0.14 and r - i=0.41±0.05 ( [Shi 
|masaku et al.|2(joT l. The fact that the red excess occurs at a 
distance from the centre that increases with colour (largest 
in sample C) may indicate that the disk is becoming more 
dominant in the redder galaxies. 

Fig. [Tl] also shows the corresponding data for the full 
sample. The mean r — i is somewhat lower but the red excess 
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Figure 10. Colour profiles along the major axis of the thin disk 
of sample T. The opening angle is 30°. A flattening of b/a=0.24 
is assumed. The bars indicate mean errors. 



is even more apparent. In the separate A-C samples the 
halo colours are too noisy to use for a rigid analysis but 
in the full sample the red excess emerges as a prominent 
structure (Fig |12| ) . In the diagram we have also plotted the 
colour profile before we correct for the effects of the PSF. We 
can see that although the correction for the PSF is highly 
significant, the red excess persists after the correction. As 
an additional support of the existence of the red excess we 
also look at the profile based on the deconvolved images. As 
we noted above, we do not trust the profiles at faint surface 
brightness levels but as is seen from Fig. |14| the colour profile 
at intermediate regions agree quite well with the data just 
discussed. The colour in the halo, r — i ~0.9, is close to the 
colours derived from the model corrected data. We finally 
note that all the thick disk subsamples, as well as the full 
samples show a weakly significant drop in r — i at very large 
distances from the centre. This may indicate that there is a 
transition from the red halo to a bluer component or that 
reduction problems become serious at these extremely faint 
levels. The uncertainties are however too large to allow for 
further analysis. 



3.4 Further tests for bias in the thick disk/halo 
colours 

We have shown that the PSF will affect the colours such 
that the colour of the faint halo regions will become redder. 
Corrections for this reddening apparently are not sufficient 
to bring down the colours to values of a normal old stel- 
lar population. We have also argued that the reduction and 
stacking procedure is not likely to have caused the observed 
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Figure 12. Colour profiles along the polar direction of sample T. 
The opening angle is 60° . A spherical shape is assumed. The bars 
indicate mean errors. The dashed line shows the colour profile 
before correction for PSF effects. The mean surface brightnesses 
in the region 100-120 pixels are /i s =30.8, /i r =30.2 and /^i=29.4 
mag arcsec - 2 . 



excess. The fact that the excess occurs at a radius where 
the galaxy is about to drown in the noise of the background 
however, made us concerned about other potential problems. 



3.4-1 Contamination by background objects 

In the process of cleaning the images free from stars and 
background galaxies to finally derive a flat background with 
a low noise level, we ran into the problem of how to deal with 
faint objects superposed on the galaxy, contaminating the 
results. It is extremely difficult to distinguish between star 
clusters and/or H II regions in the target galaxy and super- 
posed background galaxies of luminosities close to the detec- 
tion limit. The faintest background objects were therefore 
not removed from the galaxy images after rotation, align- 
ment and the second sky subtraction. Thus we treated the 
regions defined by the target galaxies differently from the re- 
gions free of bright and large objects. In the first case we did 
not remove the faint small objects while in the second case 
we did, both in the automatic identification procedure and 
in the median filtering. It would have been more consistent 
if we would have chosen not to remove the field objects, but 
that would have increased the uncertainty in the zeropoint 
of the background flux leaving us with unwanted structures 
caused by distant galaxy clustering. 

Although most of the influence of the superposed back- 
ground galaxies should be removed by the median stacking 
one may ask if it is possible that the red excess is due to 
contamination of high redshift galaxies. Indeed, for exam- 
ple an elliptical galaxy, the colours rapidly become very red 
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Figure 11. Colour profiles in the polar direction of the thick disk in the subsamplcs and the total sample. The integration is performed 
along two opposite wedges with opening angle 120°. A disk flattening of b/a=0.39 is assumed. The bars indicate mean errors. Only data 
with errors below 1 mag have been plotted. The mean surface brightnesses in the region 100-110 pixels of sample T are fi g =30.9, /i r =30.2 
and fii=29.3 mag arcsec -2 . Please note the difference in radial extension between the subsamples (A, B and C) and the full sample (T). 
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in both g — r and r — i at intermediate redshifts (see e.g. 
Csabai et al.|2003 |. Since i is the brightest part of the spec- 
trum of the redshifted galaxies we would see the excess in 
r — i first while the light from the LSB target galaxy would 
still dominate the g band, so the g — r would not be much 
affected. Moreover, if most of the light comes from galaxies 
at redshifts slightly above z=l, g — r would become bluer 
while r — i would still be very red. From this perspective, 
it would seem plausible that this is the explanation of the 
excess. However, there are counter arguments. 

In the stacking procedure we used median filtering to 
obtain the final result. If we instead had done an averaging 
of the images, the possible contribution from background 
galaxies would have increased and the colours would have 
become even redder. In Fig. [13] we show the average colours 
of the full sample. Although the scatter is somewhat larger, 
the general trend is the same as in the median filtered im- 
age. What might seem worrying is that the red excess starts 
at smaller radii in the averaged image, as we would expect 
if the background objects have more influence. Since the 
colour further out is the same in the averaged and median 
filtered image within the errors, we would have to conclude 
that most of the light in the outermost region is due to 
background objects. This is in conflict with the precondi- 
tion that SExtractor would tend to miss out objects that 
are superposed on a halo with a surface brightness signif- 
icantly above the sky level but remove them at regions of 
lower surface brightness. The result would be that the red 
excess would show up at regions of surface brightness slightly 
above the background but then fall back to normal at even 
lower surface brightnesses. There is also a second argument 
that could explain why the average colours become redder 
at a smaller radius than in the median filtered image. Since 
the galaxies are not in fact entirely similar in structure and 
stellar content, we would expect that, if a red excess exists 
at faint levels, it would appear at different radii in different 
galaxies. The average would immediately be affected when 
we reach the radius where the "first" galaxy with red ex- 
cess enters, but the median would not. Therefore we would 
expect to see the upturn in r — i at a smaller radius in the 
average filtered image. 

Perhaps the strongest argument against a background 
bias emerges when we compare the halo colour profiles Fig. 
[12] and the disk colours, shown in Fig. |10| We can clearly see 
the red excess in the halo profile but not in the disk profile 
at the same surface brightness. This cannot be explained 
by problems with background contamination. Therefore we 
conclude that it is far more likely that the red excess is 
intrinsic to the LSB galaxy sample. 

3.4-2 Inhomogeneities in the halo properties 

The small difference between the average and median stacks 
indicates that the distribution of halo-to-halo properties 
must be fairly symmetric. While halo properties may well 
vary substantially from galaxy to galaxy, the red excess can- 
not be due to a small fraction of bright and abnormal halos 
within the sample of galaxies. Moreover, the light sources 
responsible for the red colours must be relatively smoothly 
distributed throughout the halos. If the halo light were con- 
centrated into a small number of randomly positioned pixels 
within each halo, the light from such sources would be taken 
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Figure 14. The halo colour profile after PSF deconvolution. Bars 
are mean errors. 



into account in the average stack, but ignored by the median 
stack. The quantitative difference between the average and 
median properties could in principle also be used to con- 
strain scenarios in which the red halo colours are attributed 
to a limited number of massive star clusters with abnormal 



stellar initial mass functions ( Zackrisson fc Flynn|2008[ ) , but 
such an analysis is outside the scope of the current paper. 



3-4-3 Improper sky subtraction 

We finally looked at the potential problems with the sky 
subtraction. Although all the objects found by SExtractor 
were flagged, there will still be a weak remnant border of 
the transition region between objects (e.g. a star and sky). 
It is therefore an obvious risk that this remnant will pro- 
duce an artificial signal in the sky. Bright stars superposed 
on the galaxy are removed by hand and will therefore proba- 
bly produce a different remnant from that which is produced 
by the automatic stellar removal. As a consequence, the sky 
around and 'under' the galaxy would look slightly different. 
After the sky removal there would therefore be a region be- 
tween the galaxy and the sky that could show up either as a 
positive or negative remnant. But as can be seen from Fig. 
[|]and[9] there is a smooth transition from the galaxy lumi- 
nosity profiles to the sky noise. Thus, these biases must be 
too weak to significantly influence the results. 



4 WHAT IS THE CAUSE OF THE RED 
EXCESS? 

Here, we explore a number of possible explanations for the 
anomalous halo colours, involving both stellar and interstel- 
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Figure 13. Colour profiles based on the average of the stacked images of the full sample. Thick disk and halo distributions. Bars are 
mean errors. 



lar effects. We will demonstrate that while a standard halo 
population is unable to account for the colours, there are no 
less than four different mechanisms (a bottom-heavy stellar 
IMF, nebular emission, extended red emission and extinc- 
tion of extragalactic background light) that in principle can 
account for the data. However, only two of these (a bottom- 
heavy IMF and extinction of extragalactic background light) 
can also account for similar red excesses detected by others 
at longer wavelengths. 



4.1 A stellar halo with a standard initial mass 
function 

Current simulations of the formation of stellar halos around 
disk galaxies predict that the halo should start d ominating 



over the disk and bulge at fiv ~ 29 mag arcsec ( Gauthier 



et al.| [2006), which corresponds to /i; 
if one assumes g — r = 0.4, r — i = 



w 28.6 mag arcsec - 
r — i = 0.2 as expected for a 
standard, low-metallicity halo population (see Fig. |15[ ). This 
is in good agreement with the surface brightness at which 
our red halo colours turn up (/Xj ~ 28-29 mag arcsec -2 ; see 
Figs.[8|and |12[ ). However, we argue that the halo colours at 
that point are inconsistent with any kind of normal stellar 
population. 

In Fig. |15[ we compare the halo colours from the T sam- 
ple (filled circle with error bars) to the predictions fr om four 



different spectral evolutionary models: PEGASE.2 (Fioc & 



Rocca-Volmerange|1999|),|Bruzual fc Charlot| ( |2003[ ), |Marigo 
et al.| ( |2008[ ) and |Li fc Han| ( |2008[ ). Here, we have assumed 
a Salpeter IMF {AN /AM oc M~ a with a = 2.35 throughout 
the mass range 0.08-120 M©) and an exponentially decay- 
ing star formation rate SFR(t) oc exp(— t/r) with r = 1 
Gyr (suitable for an early-type system). While the Salpeter 



IMF is well-known to overpredict the number of stars with 
masses below 1 Mq in many resolved stellar populations 



(which are better represented by e.g. the Kroupa 2001| IMF), 
this mainly affects the mass-to-light ratio and has a negligi- 
ble impact on the integrated colours. For models which do 
not allow exponentially decreasing star formation rates, we 
have implemented this star formation history by summing 
up the time steps (with appropriate weights) for an instanta- 
neous burst (i.e. single-age) stellar population. While a pure 
single-age population would produce slightly redder colours 
at the highest ages, they are still nowhere near the observed 
colours. Moreover, the perfectly coeval onset and quenching 
of star formation associated with this scenario seems unreal- 
istic given the huge spatial scales (~ 10 kpc) involved in the 
halo. For each model, two different evolutionary sequences 
are plotted, corresponding to the lowest and highest metal- 
licities allowed. In a plot of r — i vs. g — r, the predicted 
colour sequences of intermediate metallicities typically lie 
inbetween those shown. 



PEGASE.2 and the 



Bruzual & Chariot 



(2003) model 



represent some of the most widely used spectral synthesis 
codes on the market, and were previously used by |Zackris-| 
son et ak| ( |2006[ > to argue that a Salpeter-IMF stellar pop- 
ulation could not account for the halo colours presented by 
Z04. For comparison purposes, we include them here as well. 
The |Marigo et aL| ( |2008[ ) model represents the current state 
of the art in the treatment of thermally pulsating asymp- 
totic giant branch stars, which are very important for the 
interpretation of near-infrared colours (e.g. |MarastonS |2005| 



Tonini et al. 20091. The Li & Han (20081 model is one of 



the few models on the market that take the evolution of 
binary stars into account. While the prescriptions and as- 
sumptions going into these models differ, it is evident from 
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Figure 15. The r — i versus g — r colours of the halo derived from sample T (filled circle with error bars), compared to the predicted 
spectral evolution of stellar populations with Salpeter IMFs and exponentially declining star formation rates (SFR(t) oc exp(— t/r), r = 1 
Gyr). Also included are the halo colours derived by Z04 for stacked high surface brightness galaxies (white triangle with error bars). For 
each model, the evolution of the lowest (dashed line) and highest (solid line) metallicities have been plotted. The markers along these 
evolutionary tracks indicate ages of 1 Gyr (circle) and 14 Gyr (star). The different panels correspond to different models: a) PEGASE.2 
( |Fioc fc Rocca-Volmerange||1999} b) |BruzuaI fc CharIot| ( |2003| l c) |Marigo et ah] ( |2008[ > d) the |Li fc Han| | |2008| ) model for binary star 
evolution. The Li &: Han model sequences are limited to ages ^ 1 Gyr and therefore do not extend to as blue r — i for low ages as the 
other models. None of these models for standard stellar populations arc able to reproduce the observed halo colours. The arrow on the 
left side of each plot represents the dust reddening vector in the case of a Milky Way extinction curve (Pei 1992) and a g-band extinction 
of A(g) = 1.0 mag. Since this arrow runs almost parallel to the age vector at ages above 1 Gyr, any reddening due to dust which may 
be present in the halo cannot reconcile the observed halo colours with a standard stellar population. 



Fig.[l5]that the differences between their predictions for the 
evolution of the g — r and r — i colours are relatively modest, 
at least at high ages. Moreover, all models - regardless of 
the metallicity - fail to explain the very red r — i colours 
of the halo from sample T. Also shown in Fig. [15] are the 
halo colours derived by Z04 for stacked high surface bright- 
ness galaxies (white triangle with error bars). These halo 
colours are marginally consistent with ours, but both sets of 
measurements are clearly anomalously red compared to the 
models. Is it then possible that this discrepancy is simply 
due to some shortcoming inherent to all models? This does 
not seem likely, since the red halos are considerably red- 
der than other old stellar populations like globular clusters 
and elliptical galaxies (see Z04), which these models have no 
problem in explaining. 

The model predictions are formally valid for redshift 
z — 0, whereas the median redshift of the T sample is z ~ 



0.06. As long as the IMF is Salpeter-like, the fc-corrections 
on the observed colours are, however, very small. Using the 
|Zackrisson et al.] ( |2001[ > spectral evolutionary code, we have 
assessed the maximum fe-corrections on the g — r and r — i 
colours to be \A(g — r)\ < 0.1 mag and |A(r — i)\ < 0.03 mag 
for & Z — 0.001, r = 1 Gyr stellar population at ages higher 
than 1 Gyr. Since these corrections are smaller than the 
observational error bars on the halo colours, redshift effects 
have no impact on our comparison of models and data. 

The arrows in the the lefthand parts of Fig. |15^ -d rep- 
resent the dust reddening vectors expected from a Milky 



Way extinction curve (Pei 19921 with g-band extinction 



A{g) = 1.0 mag. Dust reddening makes the colours shift 
in a direction which is largely parallel to that of the age 
vector at ages higher than 1 Gyr in this diagram, and there- 
fore fails to reconcile the observed halo colours with normal 
stellar populations. Estimates of dust extinction in low sur- 
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Figure 16. The estimated population mass per resolution ele- 
ment in the T sample as a function of surface brightness. The dif- 
ferent line types represent pixel masses derived in filters g (solid) , 
r (dashed) and i (dash-dotted). The gray horizontal line indicates 
the approximate initial stellar population mass below which IMF 
sampling effects become important when observing an old (^ 1 
Gyr) stellar population with a standard IMF at gri wavelengths 
| |Cervino fc Lurid iana 2004). In the surface brightness range where 
the red halo colours are derived (fii F3 28-29 mag arcsec -2 ), the 
mass per pixel is only a few times 10 3 Mq. However, surface 
brightness fluctuations due to an insufficiently sampled IMF can 
be dampened within the typical SDSS seeing disk, which may 
boost the mass relevant for the comparison by up to a factor of 
ss 10 (large arrow). The conversion between current and initial 
population mass can boost the mass by another RJ 40% (small ar- 
row). The initial mass per resolution element may therefore be a 
few times 10 4 Mq, which indicates that we are slightly above the 
surface brightness limit where IMF sampling effects are expected 
to become important. 



face brightness galaxies (e.g. Bergvall et al.p999 Rahman 
et al. 20071 indicate far less extinction than this. Current 



estimates of the opacities of halos based on studies of back- 
ground objects (which are likely to suffer more extinction 
than stars sitting inside the actual halo) also suggest much 
lower optical extinction values ( |Zaritsky|l994[ | Menard et al. 
[20091 ). Hence, it seems highly unlikely that the extinction in 
the halo could be as high as A(g) = 1.0 mag. We therefore 
conclude that the observed colours of halos around low sur- 
face brightness galaxies cannot possibly be reconciled with 
any normal type of stellar population. 



4.2 Insufficient mass sampling of a standard 
stellar population 

When comparing the observed colours of galaxies to the pre- 
dictions of spectral evolutionary models, it is customary to 
assume that the IMF is well-sampled for all stars that con- 
tribute significantly to the integrated flux. This approxima- 
tion is reasonably safe as long as the mass of the stellar 
population analyzed is M > 10 4 -10 5 Mq (assuming a stan- 
dard IMF), but starts to break down as one moves to lower 
population masses, due to substantial scatter in the num- 
ber of the high-mass stars present ( |Cervino fc Valls-Gabaud| 
2003| |Cervino fc Luridiana||2004[ ). 

While the halo regions that we analyze are very large, 



the integrated luminosity of each pixel rapidly drops as a 
function of distance from the disk. This means that the stel- 
lar mass (as derived from the luminosity) contained in each 
resolution element will eventually fall below the safe limit as 
one moves outwards. At that point, one expects the observed 
colours to start deviating from the predictions of models 
which assume a well-sampled IMF, like the ones used in the 
previous section. The tell-tale signature of this would be a 
skewed flux distribution within the halo, with a few bright 
stars sprinkled across a more uniform background of faint 
low-mass ones. When the images of many such halos are 
combined, a median stack will reject the bright outliers and 
favour pixels containing the fainter and more common stars. 
Since high-mass and low-mass stars have very different spec- 
tra, a strong colour bias will be introduced. Could this then 
be the reason for the anomalous colours that we observe? 

There are a number of arguments that make this seem 
implausible. To begin with, the observed halo colours seem 
to be insensitive to whether the images have been stacked 
by the average or median pixel flux (see Figs 



11 



12 



13) 



Since a stack based on the median will reject the extreme 
outliers in these distributions, and thereby ignore the few 
pixels which happen to contain the rare high-mass, high- 
luminosity stars, the mean and median stacks are expected 
to give rise to radically different results if the halo light 
stems from a stellar population with a standard IMF. More- 
over, high-mass stars (red giant branch and asymptotic giant 
branch stars) tend to be significantly redder (in e.g. r — i) 
than stars on the upper main sequence. Hence, one would 
expect a bias of this type to render the halo anomalously 
blue, which is the opposite of what we observe. Therefore, 
the anomalously red halo colours cannot be attributed to 
insufficient mass sampling of a normal stellar population. 

While IMF sampling effects appear not to be an issue 
for the present study, one may ask at what surface bright- 
ness limit such effects are likely to become important? To 
address this, we in Fig.[l6]estimate the stellar mass expected 
in a single pixel as a function of surface brightness within the 
stacked halo from sample T. These estimates are based on 
stellar population mass-to-light ratios (in units of Mq/Lq) 
of M/L g = 2.1, M/L r = 1.9 and MjU = 1.7, suitable for 
a 10 Gyr old, Z = 0.001, r = 1 Gyr stellar population with 
a|Kroupa|(|2001[) IMF throughout the mass range 0.01-120 



M e , as predicted by the PEGASE.2 code (Fioc & Rocca- 
Volmerange|1999[ ). Each pixel (0.396 arcsec) is moreover es- 



timated to correspond to a linear distance of 120 pc for the 
T sample. In the region where the red halo colours are de- 
rived (fj,i ~ 28-29 mag arcsec -2 ), the population mass per 
pixel is only a few times 10 3 Mq. However, surface bright- 
ness fluctuations may be smeared out over the SDSS seeing 
disk (FWHM 1.4", i.e. an area covering ~ 10 pixels), and 
this can boost the relevant mass per resolution element by a 
factor up to ~ 10 (indicated by the large arrow in Fig. 16). 



When comparing to the IMF sampling effects modelled by 



Cervino & Luridiana (20041, we moreover need to convert 



the current stellar mass (as derived from the adopted mass- 
to-light ratios) to the initial stellar mass. This can amount 
to an additional boost of up to 40% (indicated by the small 
arrow) . 

In conclusion, we may be probing a stellar population 
with an initial mass of a few times 10 4 Mq per resolution ele- 



ment. The analysis by Cervino & Luridiana (2004) indicates 
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that an initial stellar mass of > 2 x 10 Mq would be re- 
quired (for an old (age ^ 1 Gyr) stellar population observed 
in the gri wavelength range) to ensure that IMF sampling ef- 
fects are not too serious. Hence, we are probably close to the 
surface brightness limit where IMF sampling effects are ex- 
pected to become important. This means that one should in 
general avoid using median stacks when attempting to com- 
bine images of different low-redshift objects for the purpose 
of studying regions of very low surface brightness. Alterna- 
tively, many pixels within each image could be combined to 
raise the expected stellar mass per resolution element above 
the safety limit before stacking. 



4.3 A stellar halo with a bottom-heavy initial 
mass function 

The stellar IMF has long been suspected to be universal, 
but recent observational studies based on unresolved stellar 
populations have given some support to the notion that it 



may vary both as a function of environment (Hoversten & 
Glazebrook"||2008| |Meurer et al.||2009[), and as a function of 



cosmic time 



HDave|[2008 



van Dokkum 2008 Wilkins et al 



2008). Attempts to constrain the IMF in the field popula- 



tion of the LMC using star counts (e.g. Gouliermis et al 



2006) have also yielded results seemingly inconsistent with 



Zackrisson et al 



(2006) 



the standard IMF. 

Using the PEGASE.2 model, _ 
demonstrated that a bottom-heavy IMF with slope a = 4.5 
(AN /AM oc M~ a ) provided a reasonable fit to the red halos 
of both blue compact galaxies and the stacked high-surface 
brightness disks of Z04. This implies a stellar halo popula- 
tion with a huge number of low-mass stars and a very high 
mass-to-light ratio. 

Since a significant fraction of the cosmic baryons at low 



redshifts remain unaccounted for (e.g. Prochaska & Tumlin- 
son p009t , one could envision that red halos with a bottom- 



heavy IMF may serve as reservoirs for some of them. How- 
ever, Zackrisson & Flynn ( 2008 1 argue that scenarios of this 



type are inconsistent with direct star counts in the Milky 
Way halo, unless the low-mass stars are predominantly lo- 
cated in star clusters. A smooth halo population obeying this 
extreme IMF could remain viable if the surface brightness 
of the Z04 halo had been overestimated (for instance due 
to scattered light; de Jong||2008[ ), but the mass of this halo 
would then be too low to have any bearing on the missing 
baryons. 

As shown in Fig. |15[ the halo we detect in sample T is 
about 0.2 mag redder in r — i than the Z04 halo. This places 
the halo in a region of the r — i vs. g — r diagram which 
is somewhat difficult to reconcile with current models based 
on a bottom-heavy IMF of the type proposed by Zackrisson 
|et al.| ( [2006[ ). The problem is illustrated in Fig.|17| where we 
compare the halo colours observed by Z0 4 to the a = 4.5, 
Z = 0.008 PEGASE.2 model favoured by |Zackrisson et al. 



(2006). This combination of IMF and metallicity is unable 
to explain the high r — i colour of the halo from sample T at 
any age. The fit can be improved by adopting an even more 
extreme IMF (a > 4.5), but this also pushes the model track 
towards redder g—r colours. This is shown in Fig. |17b , where 
stellar populations with Z = 0.008 and a = 4.6, 4.7, 4.8, 4.9 
and 5.0 (from left to right) are shown. Even though the 
model tracks of the a = 4.8-5.0 IMFs (the three rightmost 



evolutionary sequences) are formally within the error bars 
of the halo from sample T, this fit would imply a curiously 
young halo ( < 1 Gyr old). Moreover, the required metallicity 
of Z = 0.008 is unexpectedly high compared to the stellar 
halos of the Milky Way (Z < 0.001; e.g. |Carollo et al.|2007 l 
and Andromeda (Z « 0.005 for the high surface brightness 
debris and Z » 0.001 for the smooth underlying halo; e.g. 



Chapman et al. 2006, Kalirai et al. 2006, Richardson et al. 



2009) 



Just like in the case of Salpeter-like IMFs, this compar- 
ison between data and models neglects the non-zero redshift 
of the galaxies in the T sample (median redshift z w 0.06). 
Using the Zackrisson et al. (20011 spectral synthesis code 



(which gives results similar to those of PEGASE.2 in the 
relevant r — i vs. g — r diagrams), we have assessed the likely 
fc-corrections for bottom-heavy IMF populations with pa- 
rameters identical to those used in Fig. [lTj Even though the 
fc-corrections become somewhat larger than for stellar popu- 
lations obeying Salpeter-like IMFs, they still remain smaller 
than the observational error bars on the halo colours. Hence, 
redshift corrections have no major impact on the outcome 
of our comparison. 

While current models for bottom-heavy IMFs seem 
hard-pressed when trying to fit the observed halo colours, 
the prescriptions for the evolution and spectra of low-mass 
stars used may still be too crude for the best-fit parameter 
values (Z, age and to some extent the IMF slope a) to be 
taken at face value. These stars (with masses M < 0.8 Mq) 
give a negligible contribution to the integrated flux at optical 
wavelengths in the case of a normal IMF, but become impor- 
tant when bottom-heavy IMFs are considered. As demon- 



strated by Casagrande et al. (2008), the best models avail- 



able for low-mass stars still produce temperatures (colours) 
that are considerably hotter (bluer) than observed for low- 
mass stars. While bottom-heavy IMFs are generically pre- 
dicted to shift r — i in the right direction to explain the 
red halo colours, attempting to attain an exact match to 
g — r and r — i may therefore be a futile exercise, given the 
limitations of current population synthesis codes. Hence, we 
conclude that a stellar population with an unusually high 
fraction of low-mass stars still represents a plausible ex- 
planation for colours of the halo, although this population 
probably contains no more than a tiny fraction of the miss- 
ing baryons. A number of physical mechanisms that might 
explain why the IMF appears bottom-heavy in low-density 



environments are discussed by Elmegreen (2004 1 



4.4 Extended red emission 

In principle, the red halo colours could be due to the dust 
luminescence phenomenon known as extended red emission 
(ERE). While the dust component responsible for the ERE 
has yet to be identified, ERE has been reported in a wide 
range of environments where far-ultraviolet radiation and 
dust are known to coexist, including reflection nebulae, HII 
regions, the halo region of the starburst galaxy M82 and 
the high-latitude diffuse interstellar medium and cirrus in 
the Milky Way (for a review, see |Witt fc Vijh||2004| . One 
may speculate that similar conditions could be met in the 
halos of LSB galaxies, provided that some dust is present in 
these regions and that far-ultraviolet radiation from the disk 
stars can escape in the polar direction (or, alternatively, that 
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Figure 17. The r — i versus g — r colours of the halo derived from sample T (filled circle with error bars), compared to the pre dicted 
spectral evolution of stellar populations with various IMFs, based on the PEGASE.2 model (Fioc & Rocca-Volmerange 1999k. Also 



included are the colour derived by Z04 for stacked high surface brightness galaxies (white triangle with error bars). For all model tracks, 
a metallicity of Z = 0.008 and an exponentially declining star formation rate (SFR(t) <x exp(— t/r) with r = 1 Gyr) has been assumed. 
The markers along the evolutionary tracks indicate ages of 1 Gyr (circle) and 14 Gyr (star). The arrow on the left side of each plot 
represents the dust reddening vector in the case of a Milky Way extinction curve (Pci 1992) and a <?-band extinction of A(g) = 1.0 mag. 
a) A stellar population with a = 4.5 (solid line). This model provides a reasonable fit to the halo colours derived by Z04, but not to 
the halo colours measured by us. b) Stellar populations with a = 4.6, 4.7, 4.8, 4.9 and 5.0 (solid lines, from left to right). Clearly, the 
a = 4.8—5.0 IMFs (the three rightmost lines) can fit the halo colours from the T sample, but at the cost of having an unexpectedly young 
(^ 1 Gyr) halo population. 



there are a few young stars sprinkled throughout the halo). 
The ERE manifests itself as a wide (~ 1000 A) emission 
feature in the wavelength range from 6100-9500 A (i.e. inside 
either the SDSS r, i or z bands), with the peak wavelength 
varying from object to object. At the current time, we have 
no robust way of ruling out the possibility that the i-band 
excess detected in the halo of the stacked disks is due to 
ERE. However, there are a number of arguments that speak 
against this hypothesis. 

To begin with, the peak wavelength of the ERE is 
known to correlate with the density of the far-ultraviolet 
radiation field ( jSmith fc Witt||2002[ ), and having the ERE 
peak in the SDSS i-band of sample T would (even after cor- 
recting for redshift effects) probably require a typical ioniza- 
tion parameter higher than that measured in the interstellar 
medium of the solar neighborhood. It seems reasonable that 
the typical far-ultraviolet radiation density should be lower 
than this at several kpc away from the disk of LSB galaxies. 
Studies of ERE from high-latitude clouds above the Milky 
Way (e.g. |Szomoru fe Guhathakurta|1998] |Witt et al.|2008 l 



suggest an ii-band ERE peak (at wavelengths intermediate 
between the SDSS r and i bands) for these objects - which 
may be similar to the hypothetical ERE-emitting clouds 
above LSB disks - although exceptions do exist. Most of 
the Galactic high-latitude clouds currently searched for ERE 
are moreover located within 300 pc of the Milky Way disk. 
Even if one assumes a radiation field as strong as that from 
the Milky Way, interstellar gas clouds at ~ 10 kpc (where 
the red halo colours are measured) would be expected to 
display ERE peaks at even shorter wavelengths. The ERE 
also tends to coexist with a very blue component of dust- 
scattered continuum radiation, whose relative contribution 
depends on the angle between the line of sight and the di- 
rection from which the ultraviolet photons are originating. 



Even if the ERE were to peak in the SDSS i-band, it is not 
clear that it would give rise to g — r and r — i colours sim- 
ilar to those observed. Differences in dust content between 
disk galaxies of high and low surface brightness may, how- 
ever, drive the prediction in either direction. Increasing the 
optical depth of dust could reduce the far-ultraviolet flux 
at any given distance from the disk, but at the same time 
potentially increase the density of the ERE carriers. 



As pointed out by Zackrisson et al. (20071, explaining 



the i-band excess as due to ERE, would require different ex- 
planations for the red excesses reported in this band and at 
longer wavelengths. For example, the published surface pho- 
tometry on NGC 5907 ( |Sackett et al.|1994| |James fc Casali 
p^[Rudy et al.|1997||bequeux et al ||1998| |1996| ) suggests 



that that the red halo of NGC 5907 displays an excess flux 
compared to the expectations for a low-metallicity halo pop- 
ulation with a standard IMF - in both 7, J and K bands. No 
ERE has yet been reported in this wavelength range. The 
red halos of blue compact galaxies have also been detected 
through an excess in the if-band (Bergvall fc Ostlin||2002 



Bergvall et al.| [2005 Zackrisson et al. 2006). Occam's razor 
would argue for a single explanation for the red halo excesses 
of all objects and in all filters in the wavelength range from 
I to K. On the other hand, this argument based on simplic- 
ity may be weakened if some of these detections turn out 
to be spurious. |Zheng et al.| ( |1999| ) has raised serious con- 
cerns about the reliability of the surface photometry data 



reported for the NGC 5907 halo (but see Zepf et al. 2000 



for a different view). Moreover, |Zackrisson et al.| ( |2009[ ) ar- 
gue that, due to recent developments in the field of near-IR 
spectral synthesis modelling, the observed colours of many 
of the "red halos" of blue compact galaxies may now have 
other, more mundane explanations. 
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Figure 18. The r — i versus g — r colours of the halo derived 
from the T sample (filled circle with error bars) and the Z04 halo 
(open triangle with error bars), compared to the predicted colours 
of a photoionized nebula (metallicity Z = 0.001, hydrogen den- 
sity n(H) = 10 — 3 cm -3 and filling factor / = 0.1) at redshifts 
from z = to z = 0.1 (gray squares). The solid line indicates the 
spectral evolution of a halo-like stellar population [Z = 0.001, a 
Salpeter IMF and an exponentially decaying star formation rate 
(SFR(t) oc exp(— t/r) with r = 1 Gyr) as predicted by the |Za-| 
|ckrisson et al.| ( |2001| ) model at z = 0.06 (the median redshift of 
sample T). Markers along this line indicate ages of 1 Gyr (circle) 
and 14 Gyr (star). At redshifts z ^ 0.04, this model nebula can in 
principle generate the very red r — i halo colours, but has difficul- 
ties in simultaneously explaining the g — r data. A superposition 
of nebular emission and a normal stellar halo population also fails 
in this respect. 



4.5 Nebular emission 

Another possibility is that the halo colours could be con- 
taminated by extra-planar gas photoionized by hot stars in 
the disk. This scenario was investigated by Zackrisso n et al.| 
(2006) and found to give r — i colours that were much too 
blue to be reconciled with the colours derived by Z04. This 
would imply that, even though nebular emission could pos- 
sibly account for some of the light detected at large pro- 
jected distances from the disk, correcting for this contribu- 
tion would just make the underlying halo redder and there- 
fore even more difficult to explain. The models presented by 



Zackrisson et al. ( 2006 I were, however, based on the approx- 
imation of zero redshift (i.e. no fc-correction on the broad- 
band fluxes) . In reality, the redshift of the Z04 sample peaks 
at z m 0.05, whereas the median redshift of our sample T 
is z rs 0.06. While these redshifts may seem small, they 
nonetheless turn out to have substantial effects on the pre- 
dicted colours of photoionized gas. 

Because of the high emission-line equivalent widths pre- 
dicted for a purely nebular spectrum, fc-corrections can be- 
come huge even for small redshifts. For the SDSS filters used 
here, the Ha emission line shifts out of the r filter and into 
i at z ~ 0.036, which means that the the r — i colour of 
photoionized gas turns very red at redshifts higher than 
this. This is demonstrated in Fig. |18| where we plot the 
redshift dependence (gray squares indicating the evolution 
from z — to z — 0.1 in steps of Az — 0.01) of the colours 
predicted for an HII region with low density (n(H) = 10 -3 



traplanar diffuse ionized gas (e.g. Ros sa fc Dettmax] 2000 1. 
This model nebula has been generated using the |Zackris-| 



son et al. (20011 spectral synthesis code, assuming a low- 
mctallicity (z7 star s = Z %m = 0.001). 10 Gyr old stellar pop- 
ulation with a Salpeter-IMF and disk-like star formation 
history (SFR(t) oc exp(— t/r), r — 6 Gyr) as the ioniza- 
tion source. The colours of the resulting nebula are seen to 
evolve dramatically as a function of redshift, eventually (at 
z ^ 0.05) reaching r — i colours even redder that those of 
the halo colours observed by Z04 (open triangle) and by us 
(filled circle). Hence, Zackrisson et al. (20061 may have been 
too hasty in dismissing nebular emission as the cause of the 
red halo colours. 

As the Ha emission line redshifts out of the r filter, g — r 
admittedly turns blue at the same time as, r — i turns red. 
Because of this, the nebular model presented in Fig. [18] is 
unable to simulataneously reproduce both halo colours at 
any of the relevant redshifts. Mixing this nebular spectrum 
with that of a normal halo-like stellar population (solid line 
in Fig. 18 1 does not allow the observed g — r and r — i colours 
to be reproduced either. However, predicting the exact op- 
tical colours of extraplanar nebular gas is a very difficult 
endeavour, and we cannot rule out the possibility that a 
model more realistic than the one used in Fig. [18] may be 
able to provide a better fit to the observations. 

The nebular model presented assumes a complete, 
ionization-bounded Stromgren sphere of constant density, 
whereas the halo observations are likely to probe only the 
outer regions of the ionized cloud surrounding the stacked 
disk. |Zackrisson et al!] ( |2006| ) attempted to model the spec- 
tra of individual sightlincs through the outskirts of HII re- 
gions, and predicted a large scatter in the predicted optical 
colours. The presence of shocks, turbulent mixing layers, 
density gradients and a multiphased medium would make 
the problem even more complex. One way to test whether 
nebular emission really is responsible would be to stack a 
sample of galaxies all located at z < 0.036. If a halo with an 
anomalously red r — i colour is detected in such a sample, 
nebular emission can be ruled out as the primary reason. 
Unfortunately, our T sample is not sufficiently large for this 
exercise. Attempts to remove all z ^ 0.036 objects from the 
sample results in a halo colour profile with too large error 
bars to reliably claim the detection of a red excess. 

While the nebular models used here and in Zackrisson 



et al. ( 2006 I have turned out to have a very strong redshift 
dependence in r — i, this is not the case for colours like 
V — K. Therefore, the |Zackrisson eTkl] pOOo) argument 
- that nebular emission would be far too blue to explain 
the anomalously red V — K colours reported around other 
types of galaxies - still holds. Postulating photoionization as 
the primary explanation for the flux excess detected in the 
i-band (Z04 and the present paper) would require a combi- 
nation of two different mechanisms to explain red halos in 
general. Hence, the same issues raised against extended red 
emission apply to nebular emission as well. 



4.6 The extinction of extragalactic background 
light 

One of the obvious obstacles when attempting to do surface 
photometry at surface brightness levels as faint as \jh ~ 29 
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Figure 19. The effect of EBL extinction in a diagram of r — i vs. g — r. The solid line indicates the evolution of a metal-poor (Z = 0.001) 
stellar population with a Salpctcr initial mass function and an exponentially decaying star formation rate (SFR <x exp — t/r with r = 1 
Gyr), based on Ma rigo et al.| ([2008} isochrones. White markers represent ages of 1 Gyr (circle) and 14 Gyr (star). The filled circle 
represents the colours of the halo detected in the T sample (at fii RJ 28-29 mag arcsec -2 ). Dashed lines indicate how a spatially constant 
EBL extinction would shift the observed colours away from the model predictions under the assumption of a Milky Way extinction curve. 
The differently-sized triangles along these tracks correspond to surface brightness levels of in = 27, 28 and 29 mag arcsec -2 (largest to 
smallest triangles), a) Results assuming an EBL extinction of A(g) = 0.15 mag and fiducial diffuse EBL levels Zackrisson et al. (2009). At 
in = 29 mag arcsec -2 , EBL extinction can indeed bring the colours of a normal stellar population into agreement with those observed. 
This particular solution does, however, require a fairly young (« 3 Gyr) stellar halo population, b) Results assuming an EBL extinction 
of A(g) = 0.03 mag and diffuse EBL levels slightly offset from the fiducial ones (see main text for details). In this case, EBL extinction 
effects can reconcile a 10 Gyr old stellar halo with the colours observed at m 28-29 mag arcsec -2 . 



mag arcsec -2 is the challenge of subtracting the sky with suf- 
ficient accuracy. The night sky has a typical surface bright- 
ness of ^ ~ 20.3 mag arcsec -2 at the site of the SDSS 
telescope (Apache Point), which means that the signal we 
are measuring is about 3000 times fainter than the sky. To 
measure a surface brightness with an error of ±0.1 mag at 
such isophotal levels requires that the sky can be subtracted 
with an accuracy of about one part in 30000. While we have 
argued in section 3 that the night sky can indeed be mea- 
sured with the required accuracy by stacking ~ 1500 SDSS 
images, there may be a subtle systematic problem involved 
in the actual subtraction which we cannot fully correct for. 



This problem (first recognized by Zackrisson et al.|2009[ ) 
stems from to the fact that while most of the sky flux orig- 
inates from regions between the telescope and the galaxies 
that we target, a small fraction - the extragalactic back- 
ground light (EBL) - comes from behind. The EBL at opti- 
cal wavelengths is believed to be the product of direct and 
reprocessed starlight emitted over the entire star formation 
history of the Universe, and hence stems from objects at 
vastly different redshifts. Provided that the existing EBL 
measurements are correct, most of this light appears to be 
diffuse (i.e. unresolved with all existing instruments). Unlike 
the other components of the night sky flux, the EBL can be 
subject to extinction by dust present in the halos of our 
target galaxies, thereby invalidating an implicit assumption 
adopted in all current surface photometry measurements, 
namely that the sky flux and the flux from the target ob- 
jects are unrelated. The extinction of the EBL may lead to 
a slight depression of the overall sky flux across the galaxies 
(and their halos) in our sample. Since our sky subtraction 
method is based on estimates of the sky flux well away from 



the target galaxies - where the EBL suffers less extinction 
than in the regions where the halo colours are measured - 
a systematic oversubtraction of the sky may result. As ar- 



gued by Zackrisson et al. (20091, even very small amounts 



of extinction can in principle lead to spurious features in ra- 
dial surface brightness profiles and colour maps of extended 
objects. 

|Zackrisson et ak] ( |2009[ > showed that even though the 
halo colours derived by Z04 could in principle be attributed 
to EBL extinction effects, this requires both an unresolved 
EBL level at the high end of what is allowed by current 
error bars, and an EBL extinction of A(g) = 0.2 mag in the 
region where the halo colours are measured (w 8 kpc above 
the disk in the Z04 sample). Our halo colours are, however, 
derived at fii ~ 28-29 mag arcsec -2 , which is about 1-2 mag 
fainter than that of the region where Z04 note the onset of 
anomalous colours (/ij w 26.7 mag arcsec -2 ). This means 
that a similar colour shift can be produced for a somewhat 
lower EBL level and/or a lower halo opacity. 

This is demonstrated in Fig.|19| where we have adopted 



the formalism from Zackrisson et al. (2009) to predict the 



colour shifts in a diagram of r—i vs. g—r resulting from a fail- 
ure to account for EBL extinction. The solid line represents 
the predicted evolution of an old, metal-poor (Z = 0.001) 
stellar halo population modelled using |Marigo et aT] p008| 
isochrones with a Salpeter initial mass function. The dashed 
lines in Fig. |19^ indicate how the observed colours would 
shift due to EBL extinction if we assume an extinction of 
A(g) = 0.15 mag in the halo. The shift becomes increasingly 
severe at fainter surface brightness levels (triangles along the 
dashed lines indicating the colours at the — 27, 28 and 29 
mag arcsec -2 isophotes). As discussed by Zackrisson et al., 
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the direction of the shift depends on the intrinsic colours of 
the halo, which is why three different dashed lines are in- 
cluded. Indeed, one of the dashed lines passes very close to 
the observed colours at /ii w 29 mag arcsec -2 . No artificial 
boost of the observed EBL level (like that used by Zackris- 
son et al. to fit the Z04 data) is required. This particular line 
originates from a fairly young (« 3 Gyr) halo population, 
but similar fits for older populations can be achieved by as- 
suming slightly different colours for the EBL (well within the 
observational errorbars). This is shown in Fig. 19 3, where 
we have adopted fi g = 26.1 mag arcsec -2 , fi r = 25.1 mag 
arcsec -2 and f^i = 25.4 mag arcsec -2 for the diffuse EBL 
(0.2 mag fainter in gi and 0.3 mag brighter in r than the 
fiducial values). In this case, the observed colours are re- 
produced using a 10 Gyr old halo population and a halo 
extinction of just A(g) = 0.03 mag. 

The anomalous colours of the halo in sample T turn 
up at « 80-110 pixels from the plane of the disk (see 
Fig. 12 1, which converts into w 10-13 kpc (assuming 120 



pc per pixel). It is difficult to judge whether an extinction 
of A(g) = 0.03-0.15 mag is reasonable in this region, since 
no studies of the opacity profiles of halos have specifically 
targeted low surface brightness galaxies. However, |Zaritsky] 
( |1994[ ) and |Menard et al.| ( |2009| ) have measured optical ex- 
tinctions in the « 0.03-0.15 mag in the halos of other types 
of galaxies (at even greater projected distances than what 
we probe here). Dust is also well known to exist several kpc 



above the plane of edge-on discs in general (e.g. Alton et al 



2000) 



The dust column density at the site of our halo obser- 
vations can be estimated from 



'S'dust — 



0AA 



K ext log(e) 



(1) 



where A is the extinction in magnitudes, K cxt is the dust 
extinction coefficient and Sdust is the dust column den- 
sity. The dust extinction coefficients in the V band for the 
Milky Way, L MC and SMC are given by |Weingartner fc| 
Draine (2001 1. Since the effective wavelength of the g band 
is close to that of Johnson B, we will assume E(g — V) = 
E(B - V) ~ A(V)/3. Thus A{g) = 4A(V)/3. We expect 
the metallicity in the halo to be low and therefore the most 
appropriate value to use here is that of SMC, K cxt (V) = 
1.54 x 10 4 cm 2 g -1 and thus K e *t(g) = 2.05 x 10 4 cm 2 g -1 . 
A(g) = 0.03-0.15 mag then corresponds to dust column den- 
sities of 1.4xl0 -6 - 6.7xl0" 6 g cm . Assuming a pixel 
size of 120 pc, these numbers correspond to a a range in 
dust mass surface density of Ed us t = 90-460 Mq pixel - . 
We may compare these values with the stellar mass den- 
sity per pixel that we derived in Section 4.2. In the g band, 
at a surface brightness of fi ~ 29 mag arcsec -2 , we ob- 
tain E stars ss 3.6 x 10 3 Mq pixel -1 . Thus the dust-to-stellar 
mass density ratio is Md us t /Mstars ~ 2-13%. How does 
this compare to the global dust properties of LSB galax- 
ies? Rahman et al. (20071 recently presented dust masses 



and total stellar luminosities for a few LSB galaxies. From 
these data, we estimate a typical dust-to-stellar mass ratio 
of Mdust /Altars ~ 1% by assuming a stellar mass-to-light 
ratio of M/L » 1 ( |Zackrisson et al.|2005 l. 

This indicates, that to accomodate a halo extinction of 
A(g) = 0.03-0.15 mag, the opacity would have to decline 
slower than the stellar surface brightness as one moves in 



the polar direction away from the disk. A dust distribu- 
tion more extended than that of stars has already been con- 



firmed within high surface brightness disks (e.g. Holwerda 
et al.||2005| |Munoz-Mateos et al.||2009|). Considering that 



halo opacity profiles (Zaritsky 1994; Menar d" et al.|2009 l ap- 



pear to fall off much slower than halo surface brightness pro- 
files (Z04 and this paper) would suggest a similar trend for 
the dust-to-stellar mass ratios in halos, although the selec- 
tion criteria of these different studies do not exactly match. 

If the red halo colours are due to EBL extinction ef- 
fects, one may naively expect the thin disk to display a sim- 
ilar colour anomaly at comparable surface brightness levels. 
However, the colour profiles along the thin disk (Fig. \10\ in 
sample T show no obvious sign of this, and neither do the 



disk profiles presented by Z04. Zackrisson et al. (20091 dis- 



cuss two possible explanations. Firstly, the colour shifts are 
very sensitive to the intrinsic colours of the affected stellar 



population (as seen in Fig. 19 1), and a disk population is not 
expected to have colours identical to that of a stellar halo. 
Moreover, the extinction curve may be grayer in the disk 
(e.g. Holwerda et al.||2005 |. Both of these effects may con- 
tribute to diminish or alter the offset vector of any colour 
shifts related to EBL extinction. Hence, the lack of a red 
excess in the disk does not rule out EBL extinction as the 
cause of the red excess in the halo. As argued by |Zackrisson| 
et al. ( 2009[), EBL extinction may also be able to explain 



the anomalous halo colours detected around other types of 
galaxies and at longer wavelengths. 



5 SUMMARY AND CONCLUSIONS 

We present the results of an investigation of the statistical 
photometric properties of a homogeneous set of 1510 LSB 
disk galaxies drawn from the SDSS DR5. There have been 
claims that the faint outskirts of normal disk galaxies and 
blue compact galaxies tend to show colours that are redder 
than predicted by stellar evolutionary models - the so called 
'red halo' phenomenon. We have therefore focused on deriv- 
ing the characteristic colours of the faintest regions in the 
halos of our LSB galaxies. For this purpose we use a sta- 
tistical approach to reach the faintest levels. Thus we have 
rescaled, aligned and stacked the g, r and i images, thereby 
increasing the S/N of the combined image. In the finally 
stacked images we typically reach a surface brightness of 
fx ~ 31 mag arcsec -2 . 

We derive surface brightness and colour profiles based 
on median filtering and averaging of the images in the stack. 
When we reach fj, ~ 28-29 mag arcsec -2 , a significant red 
excess appears in r — i. We find no support for the notion 
that this effect can be attributed to scattered light. After 
a careful assessment of PSF effects, the combined colours, 
g — r — 0.60 ± 0.15 and r — i — 0.80 ± 0.15, are incompatible 
with colours of a normal stellar population. The fact that no 
similar colour anomaly is seen at comparable surface bright- 
ness levels along the disk rules out a sky subtraction residual 
as the source of the extreme halo colours. We discuss possi- 
ble alternatives to explain the red colours. These include a 
bottom heavy IMF, nebular emission, extended red emission 
and dust extinction of extragalactic background light. We 
find two of these scenarios to be broadly consistent with our 
results and other claims of a red excess in the optical/near- 
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IR. The first is extinction of extragalactic background fight, 
which wouid cause a cofour-dependent oversubtraction of 
the sky background in the region where dust is present. In 
this case, a hafo opacity of A(g) = 0.03-0.15 mag at a sur- 
face brightness of fi g ~ 29 mag arcsec -2 would be needed 
to explain the observations. The second alternative is an 
abnormal IMF, dominated by low-mass stars. 
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